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Abstract 

We make use of the complex implicit representation in order to provide a 
deterministic algorithm for checking whether or not two implicit algebraic 
curves are related by a similarity, a central question in Pattern Recognition 
and Gomputer Vision. The algorithm has been implemented in the computer 
algebra system Maple 2015. Examples and evidence of the good practical 
performance of the algorithm are given. 


1. Introduction 

A central problem in Pattern Recognition is to classify a given object. If 
the object to be studied is represented as a curve, the problem reduces to 
comparing this curve with other curves, previously classihed and stored in a 
database. However, the position and size of the curves in the database do not 
necessarily coincide with the position and size of the curve to be classified. 
Therefore, from a geometric point of view we need to check whether or not 
there exists some similarity, i.e. the composition of a rigid motion and a 
scaling, transforming the given curve into some curve in the database. In 
Gomputer Vision we hnd a close problem, called pose estimation, which has 
to do with identifying a certain object in a scene. Geometrically, the problem 
is equivalent to recognising a same object in two different positions. In turn. 
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this amounts to checking if two objects are the same up to the composition 
of a translation and a rotation, so that the scaling factor does not change. 

These questions and related problems have been attacked in the applied 
mathematics literature using very different approaches: B-splines |8], Fourier 
descriptors [Hj, complex representations (TB], Moebius transformations in 
the parameter space [T], statistics [7l[ini[I2] (also in the 3D case), moments 
[El [IE], geometric invariants [SniESlEl], Newton-Puiseux parametrizations 
|13j . or differential invariants [a El Eg. The interested reader may consult 
the bibliographies in these papers to hnd other references on the matter; the 
list is really very long. 

In the above references the input is typically a “fuzzy” image, possibly 
noisy, with occluded parts, etc. and quite often a point cloud. A com¬ 
mon strategy is to adjust hrst an algebraic curve to the input, so that the 
comparison is carried out between algebraic curves. Although in some cases 
(B-splines [8], Fourier descriptors [H]) the curves are described by parametric 
representations, in most approaches the curves are represented by means of 
implicit algebraic equations. Therefore, the question of detecting similarity 
between implicit algebraic curves comes up. 

Since in the given references the input is usually noisy and the repre¬ 
sentation of the object as an algebraic curve is only an approximation, the 
comparison made between the curves is also approximate. In particular, sta¬ 
tistical and numerical methods are massively employed. However, in this 
paper we address the problem from a different perspective. We assume that 
the curves are given in exact arithmetic, and we provide a deterministic, not 
approximate, answer to the question whether or not these two curves are 
similar. In the affirmative case, our algorithm can also hnd the similarities 
transforming one into the other. The motivation for this point of view comes 
from computer algebra systems. Assume that a database of classical curves 
is stored in your favourite computer algebra system. Using the algorithms in 
this paper, the system can recognise a certain curve introduced by the user 
as one of the curves in the database. This way, the user can identify a curve 
as, say, a cardioid, a lemniscata, an epitrochoid, a deltoid, etc. 

In order to do this, we use the complex representations of the curves 
to be compared. Complex representations have already been used in [161 
[T9] , where the pose-estimation problem is addressed, and in [9], where the 
computation of the symmetries of an algebraic planar curve is studied. In 
[TB] the complex representation is exploited and combined with numerical 
strategies in order to give an approximate solution to the problem. In [19] 
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the transformation to be sought is decomposed into a rotation, which is 
approached in a deterministic way, and a translation, which is approximated. 
The approach in |9] is deterministic, but they treat a different problem. 
In our case, we use the complex representation of the curves to translate 
the problem into a bivariate polynomial system of equations with complex 
coefficients. Additionally, except in certain special cases we can also add 
a univariate equation to the system (corresponding to the rotation angle), 
which allows to solve the system in a fast and efficient way. 

Besides complex representations, we also employ usual tools in the held 
of computer algebra, like resultants, gcds and polynomial system solving. 
The resulting algorithm has been implemented and tested in the computer 
algebra system Maple 2015. Our experiments conhrm that the algorithm is 
efficient and fast. 

The structure of the paper is the following. Preliminaries are introduced 
in Section |2l In Section [3] we discuss the two cases that must be distin¬ 
guished. The most general case is addressed in Section and a special 
case is studied in Section |5l Both in Section lU and Section |5] we focus on 
orientation-preserving similarities. Although the orientation-reversing case is 
analogous, some observations for this case are provided in Section]^ Finally, 
examples and experimentation details are provided in Section 


2. Preliminaries 

Throughout the paper, we consider two plane, real, algebraic curves Ci,C 2 C 
implicitly dehned by two real, irreducible polynomials f{x,y),g{x,y) of 
the same degree n with rational coefficients. Furthermore, we will write 


y) = fn{x,y) + I/) + ■ ■ ■ + /o(x, y). 


where for p 
P of /, 


Similarly, 


0,1,..., n, fp{x, y) represents the homogeneous form of degree 

p 

fp{.x,y) = '^ttp-jjxP-^yP 

j=0 


p 

9ix,y) = gn{x,y) +gn-iix,y) H- \-go{x,y), gp{x,y) = ^bp-jjxP~^yP 

j=0 
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We will assume that Ci,C 2 C are real curves, i.e. that they contain 
inhnitely many real points, and that they are not either lines or circles. We 
have found it useful to identify the Euclidean plane with the complex plane. 
Through this correspondence, (x,y) z = x + iy and 

z + z z — z _ , 


where z denotes the conjugate of the complex number z. By performing the 
change Q into the implicit equation f{x,y), we reach a polynomial F{z,z) 
that can be regarded as a “complex” implicit equation of Ci, since it is, as it 
also happens with the usual implicit equation, unique up to multiplication by 
(possibly complex) constants. We will refer to F{z, z) as the complex implicit 
representation of Ci. The polynomial F(z, z) can be written as 

F{z, z) = Fn{z, z) + Fn_i{z, z) -h Fq{z, z), 

where Fp{z, z), p = 0,1,... ,n, is an homogeneous polynomial of degree p in 
the variables z, z, i.e. for p > 1, 

p 

Fp{z,z) = '^ap_jjzP-^z^. 

j=0 

Furthermore, one can check that for p = 0,1,...,n, fp{x,y) gives rise to 
Fp{z, z) and conversely. We will say that a term in F(z, z) has bidegree 
{n — k, k) if it can be written as ^ ■ z'^~^z^, with G C. In a similar way, by 
applying 0 into g{x,y) we reach G{z,z), where 

p 

G{z,z) = Gn{z,z) + Gn-iiz,z) - ^Gq{z,z), Gp{z,z) = '^/3p_jjzP~^Z^. 

j=0 

An (affine) similarity of the plane is a linear affine map from the plane to 
itself that preserves ratios of distances. More precisely, a map h : —?• 

is a similarity if h{x) = Ax + b for an invertible matrix A G and a 

vector 6 G and there exists R > 0 such that 

\\h{x) - h{y )\\2 = R- ||a; - yh, x,y e ]R^ 

where 11-112 denotes the Euclidean norm. We refer to R as the ratio of the 
similarity. Notice that if R = 1 then h is an (affine) isometry, i.e. it preserves 
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distances. Planar isometries are completely classified [S]. Under the compo¬ 
sition operation, similarities form a gronp, and isometries form a snbgronp 
inside of this gronp. Furthermore, any similarity can be decomposed into an 
isometry and a uniform scaling by its ratio R. We say that Ci, C 2 are similar 
if and only if there exists a similarity h such that h{Ci) = € 2 - Notice that if 
Ci,C 2 are similar, since by hypothesis they do not have multiple components 
their degrees must coincide. This explains why we assumed the equality of 
the degrees of Ci,C 2 at the beginning of the paper. Furthermore, whenever 
Cl, C2 are not lines or circles the number of similarities between them is hnite 
[T]; in fact, if Ci, C 2 are not symmetric then there exists at most one similarity 
mapping one onto the other [T]. 

A similarity can either preserve or reverse the orientation. In the hrst case 
we say that it is orientation-preserving, and in the second case we say that it 
is orientation-reversing. By identifying and C, any orientation-preserving 
similarity can be written as h{z) = az -\- b, where a,b E C. Furthermore, 
any orientation-reversing similarity can be written as h{z) = az -\-b. In each 
case, its ratio R = |a| 7 ^ 0. The next result shows the relationship between 
the complex implicit equations of two similar curves Ci,C 2 . 

Theorem 1. Let Ci,C 2 be two algebraic curves with complex implicit egua- 
tions F{z,z) and G(z,z) with the same degree. The curves Ci,C 2 are related 
by an orientation-preserving similarity h{z) = az b, if and only if there 
exists A G C, A 7 ^ 0 such that 

G{az b,az -\- b) = X ■ F{z, z) (2) 

for all z E C. Also, Ci,C 2 are related by an orientation-reversing similarity 
h{z) = az + b, if and only if there exists A such that 

G{az b,az b) = X ■ F{z, z) (3) 


for all z E C. 

Proof. We prove the result for the orientation-preserving case. Similarly for 
the orientation-reversing case. Here we identify each point (x, y) E with 
the complex number z = x -E iy. (^) If Ci, C 2 are related by an orientation¬ 
preserving similarity h{z) = az + b, then for every point z E Ci the point 
oj = h[z) belongs to C 2 , i.e. G{uj,uj) = 0. Since oj = az + b the curves 
dehned by F[z, z) = 0 and G{az -\-b,az + b) =0 have inhnitely many points 


5 


in common. By Bezout’s Theorem and since they have the same degree, ([^ 
follows. (<^=) If Q holds then whenever 2 ; G Ci, i.e. whenever F{z,z) = 0, 
then G{az + b,az + b) = 0, i.e. h{z) is a point of € 2 - Therefore, the image 
h{Ci), which is an algebraic curve, must be a component of € 2 - Since C 2 and 
h{Ci) have the same degree C 2 = h{Ci). □ 

Notice that in the orientation-preserving case, h{z) = az + b is the com¬ 
position of: (1) a rotation around the origin, of angle equal to 9, where 
a = |a| ■ e*®; (2) a scaling, where the scaling factor is |a| and the center is 
the origin; (3) a translation of vector b. In the orientation-reversing case, 
h{z) = az + b is the composition of: ( 1 ) a symmetry with respect to the x- 
axis; ( 2 ) a rotation around the origin, of angle equal to 9, where a = |a| ■ e*®; 
(3) a scaling, where the scaling factor is |a| and the center is the origin; (4) 
a translation of vector b. 

3. Similarity detection: discussion of cases. 

The equations (|^ and (|^, jointly with the conditions a, A 7 ^ 0 are nec¬ 
essary and sufficient for Ci,C 2 to be similar. Equaling the terms in z'^~^z^, 
0 < j < m, 0 < m < n at both sides of (|^ leads to a polynomial system 
S of at most N = equations with complex coefficients in the variables 
a, a, b, b, A. The solutions to that system with a, A 7 ^ 0 correspond to the 
similarities between Ci,C 2 preserving orientation. Similarly, by considering 
(|^ we reach another polynomial system S' whose solutions correspond to 
the similarities reversing orientation. Thus, Ci,C 2 are similar if and only if 
either S or S' is consistent and has a solution with a, A 7 ^ 0. In the sequel, 
we will show how to efficiently check the consistency of S. For S' the analysis 
is similar; we will make some observations on this case in Section 

Since Ci has degree n, there exists j G {0,1,..., n} such that an-jj 7 ^ 0. 
Let us see that we can always hnd j ^ n with this property. We need the 
following lemma. This result is mentioned without a proof in |T 6 ] and [9] , so 
for completeness we include a proof here. 

Lemma 2. For j G {0,1,..., n} we have exn-jj = aj^n-j, Pn-jj = Pj,n-j- 

Proof. We prove the result for the a’s; the proof for the /3’s is analogous. 
Performing the change ([^ in each term an-k,kx''^~^y^■, where /c G {0,1,..., n}. 
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we get 


iz+z\^~^ (z—z\^ 

an-k,k (^) ■ (^) = 

If we expand this multiplication, whenever ji + j 2 = j we obtain a term in 
therefore, 


^n — k,k 


s;^n-k /n-k\ n-k-ji zh 


Oin-j,j 


n 


E 


^n—k,k 
2n . j^k 


E 

ji + j2 = 3 
ji<n- k, j 2 < k 


fn — k\ 

V Ji J 



( 4 ) 


Similarly, in order to hnd ctj^n-j we substitute ([^ into an-k,kx"' ^y^-, but we 
proceed in a slightly different way: 

(z±z\'^-^ ( -z+z \^ _ 

^n-k,k \ 2 J ' \ 2i ) 

T.i.0 ■ 

When ii + £2 = j we get hence. 


_ an-k,k '^n-k /n-k\ -n-k-ii h 
2n.jfc ■ Z^4=0 V £i ^ 


C^j,n-j 


n 


E 


^n—k,k ' ( 1 ) 

2n . j^k 


E 

£ 1 + £ 2 = j 
£i < n — k ,£2 < k 



( 5 ) 

Finally, comparing the terms in Q and ([^ we notice that when k is even 
we get the same real number, and when k is odd we get two opposite purely 
imaginary numbers. So the result follows. □ 


If Q;o,n 7 ^ 0 then from Lemma we have that an,o 7 ^ 0 too; therefore there 
always exists j G {0,1,..., n — 1} such that an-jj 7 ^ 0. Now if an-jj 7 ^ 0 and 
f^n-jj = 0, the equality (|^ cannot be satisfied, and therefore Ci, C 2 are not 
related by an orientation-preserving similarity. Thus, in the sequel whenever 
o.n-j,j 7 ^ 0 we will assume fin-jj 7 ^ 0 too. Since 

G{az b,az + b) = ■ ■ ■ + f3n-j,j{ciz + b)"‘~^{az + by + ■ ■ ■ , 
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( 6 ) 


by equaling the coefficients of z"- ^at both sides of (|^, we get 

^ _ /3n-j,ja"'~^a^ 

C^n-j,j 

Now let us compare the coefficients of z'^~^~^z^ at both sides of ([^. In order 
to hnd the coefficient of z'^~^~^z^ in G{az + b,az + b), we observe that the 
terms of G{z, z) contributing to the term z"‘~^~^z^ in G{az + b,az + b) are 
the terms with bidegrees {n — j,j), — j — 1, j + ^) and (n — j — 1, j). After 

substituting z ^ az + b in these terms, the coefficient of z'^~^~^z^ at the left 
hand-side of (|^ is 

(n - + (j -f- l)Pn-j-ij+i ■ a^~^~^a^b + (3n-j-i,j ■ 


Notice that this expression does not make sense when j = n, which is the 
reason why we needed j G {0,... ,n — 1}. The above coefficient has to be 
equal to the coefficient of z"‘~^~^z^ in X-F{z, z), namely X-On-j-ij- Therefore, 
after taking ([^ into account and dividing by we deduce that 


(n j)Pn-j,j ■ h + (j + l)Pn-j-l,j+l ■ b -\- j3n-j-lj 




^n-jj 


a, (7) 


which is linear in a,b,b. By conjugating Q, we also deduce the following 
relationship: 


(j + l)Pn-j-l,j+l -b+in- j)Pn-j,j ■ b + 






a. (8) 


Let Aj = {n - j)^ ■ \/3n-jj\^ - {j + 1)^ • |/3n-i-i,j+iP- The following 
observation is crucial for us: if Aj ^ 0 , then we can write 6 as a linear 
expression b = T](a,a) from the linear system consisting of Q and (|^. 
Hence, we will distinguish the following cases. Here we assume that cun-jj ^ 
0 iff Pn-jj 7 ^ 0 , since otherwise we know that Ci and C 2 cannot be similar. 


• We will say that we are in the general case for C 2 , if there exists some 
j G {0,..., n — 1} with an-j,j 7 ^ 0 such that Aj 7 ^ 0. 

• We will say that we are in the special case for C 2 if for all j G {0,..., n — 
1 } with an-j,j 7 ^ 0 , we have Aj = 0 . 





The general or special cases for Ci can be introduced similarly. 

In order to quickly detect whether or not we are in the general or the 
special case, we can proceed as follows. Let j be the minimum of the j G 
{0,..., n — 1} such that ctn-jj, Pn-jj 7^ 0 (observe that j < [fD. Then the 
following result holds. 

Lemma 3. Assume that cxn-j,j ^ iff fn-j,j 7^ 0. The special case for C 2 
occurs if and only if \(dn-j,j I = ' \fn-j,j \ for j = 0, ... ,n — l. In particular, 

if 3 f ^ then the special case for C 2 does not occur. 

Proof. («y=) Since = (”) ■ \(dn-j,j\ for j = 0,..., n - 1, we have 



Since 


0:0 


''^6 can easily see that Aj 


dehnition (dn-jj 7 ^ 0, we have fdn-j,j 7 ^ 0 for j 
(y.n-j,j 7 ^ 0 for j = 0 ,..., n — 1 . 


0. Furthemore, since by 
0,..., n — 1. Therefore, 


(=>) By induction on k, one can prove that \l3n-{j+k),j+k\ = (j+fc)' l/^n-jyl 
for 0 < A; < n—j. Thus, when k = n—j we deduce that |/3o,n| 7 ^ 0- Therefore, 
/do,n 7 ^ 0 and by Lemma (dnfl 7 ^ 0. Hence j = 0. Furthermore, since by 
dehnition (dn-jj 7 ^ 0, we also deduce that all the an-j,fs are nonzero. □ 


Corollary 4. IfCi,C 2 are similar, the special case for C 2 occurs iff the special 
case for Ci occurs too. 

Proof. If Ci,C 2 are similar then from Equation (|^, we have that 

■ (dn-j,j = A ■ On-jj (9) 


for j = 0,..., n. Furthermore, by taking the absolute value in (|^ and hxing 
3 = Jy we get that 


|A| = 


\ldn 

-id 

■ 1 

a 

n 


^n— 

id 

1 


( 10 ) 


So from and (10), we get 


I ^ri-j,j \ 






Then the statement follows from Lemma [H 


□ 
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If the special case happens for one of the cnrves bnt not the other, from 
Corollary the cnrves cannot be similar. So from now on we will assnme 
that if the special case happens for one cnrve, it also happens for the other 
cnrve, and we will refer to this sitnation as the special case. In the next two 
sections, we provide algorithms to check similarity in each of these cases. 

4. The general case. 

By using the expressions (§, 0 and (|^ we can write X,b,b in terms 
of a, a. Therefore, the system S gives rise to a bivariate system in a, a 
with complex coefficients, and Ci and C 2 are similar iff this system has some 
solution with a 7 ^ 0. However, we will see that for a generic pair of curves 
Cl, C2, we can do better. The rough idea is to move from the bivariate system 
in a, a to another bivariate system in two new variables, which contains at 
least one univariate equation. The fact that the system contains a univariate 
equation makes the analysis of the system much easier and faster. This idea 
does not work in certain special, bad situations, described below. 

Let a = \a\ ■ e*® = |a| ■ (cos^ + isin 6 *), 9 G (— 7 r, 7 r]. We need to consider 
separately the cases cos( 6 ') 7 ^ 0 and cos( 6 *) = 0. We focus hrst on the case 
cos 6 * 7 ^ 0 ; the case cos 6 * = 0 is easier and will be treated at the end of the 
section. Under the assumption cos 0 7 ^ 0, we can write a = r- (1 Tio;), where 
r = \a\ ■ cos 9 and u = tg(9). Now observe that if Ci and C 2 are related 
by a similarity h(z) = az + b, then h maps the curve dehned by the form 
of maximum degree of Ci, fn{x,y), into the curve dehned by the form of 
maximum degree of C2, gn{x,y). The former is mapped onto the latter by 
the mapping h[z) = a^;, which is also a similarity. The curves fn{,x, y) = 0 
and gn{x,y) = 0 are the union of n (possibly complex) lines, counted with 
multiplicity. Let us denote the lines corresponding to the curve fn{x,y) = 0 
by Cl,..., Cn, and the lines corresponding to the curve gn{x,y) = 0 by 
M-i,... ,Mn. If Cl and C2 are similar, the Wfj’s are the result of rotating 
the Cj’s about the origin by an angle equal to 9. Let us see that in most 
situations we can either directly hnd cj, or hnd a polynomial P(t) such that 
a; is a real root of P(t). In order to hnd P(t), we distinguish two cases: 

• Case (1): gn{x,y) = x’^{ax + by)^ k ^ 0, k ^ n, b ^ 0. If the curves 
Cl and C2 are similar then fn{x, y) must have two irreducible factors as 
well, with multiplicities k,n — k. If x does not divide fn{x,y) we can 
move to the second case. Otherwise fn{x, y) = x^{cx + dyY~^, where 
k = ioTk = n — i, Now we have several subcases: 
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— If a = c = 0 then fn{x, y) = x^y"'~^ and gn{x, y) = x^y'^~^. In this 
case, we proceed as follows: (1) if /c 7 ^ n — k, then ii k = i then 
a; = 0, i.e. it is a root of P(t) = t, and if k i then cos 6 * = 0, 
which has been postponed; {2) if k = n — k, then either ca = 0 or 
cos 6 ' = 0 , and both possibilities must be explored. 

— If a 7 ^ 0 or c 7 ^ 0, then: 

* If k = i and n — kj^k, then ca = 0, i.e. ca is a root of P(t) = t. 

* If k i then u = —d/c, i.e. ca is a root of P(t) = t + d/c. 

* lfk = n — k = i then a; is a root of P(t) = t(t + d/c). 

• Case (2): case (1) does not occur. Here, both /„ and Qn have 
factors corresponding to lines y = mx. In the case of fn{x,y), the 
m’s are the (possibly complex) roots of p{y) = /„(!,?/); in the case of 
gn{x,y), the m’s are the (possibly complex) roots of q{y) = gn{f,y). 
Notice that 6 * = 72 — 71 , where 71 is the angle formed by one of the C/s 
and the positive x-axis, and 72 is the angle for the corresponding Aii. 
Since 72 = 7 i + 0, denoting m = tg( 7 i), m = tg( 72 ), from the tangent 
angle addition formula we get that: 

m + a; 

m = -. 

1 — mu 

Notice that mu 7 ^ 1, because we can assume that we are relating two 
linear factors of fn{x, y) and 5 'n(a^, y) different from x. Now m is a root 
of qiy) and m is a root of piy). Therefore, writing 

:= 9 , ( 11 ) 


a; is a root of 


Pit) = Resj^ ipiy), num(Q(t, y))) = 0, (12) 

where “num” denotes the numerator of the expression in brackets. 

Continuing with the description of the algorithm, since we have h = 
r]ia,a) with a = r- (l + ica) and rj linear, we get b = r]ir,u), b = r/(r, ca). By 
substituting these expressions into the system S, we reach a bivariate system 
in the variables r, u. By adding the equation Piu) = 0 to this system, the 
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curves Ci,C 2 are similar if and only if this new system has a solution with 
r G M, r 7 ^ 0 and a; G M. 

One can wonder if the polynomial P(t) in (12) can be identically zero, in 
which case the advantage of having a univariate equation is lost. The answer 
is affirmative, as shown by the following proposition. 


Proposition 5. The polynomial P{t) in (12) is identically zero iff p{y) and 
q{y) are both divisible by y‘^ + 1. 


Proof. Let 


Hence 


which implies 


p{y) = Yl{y (i{y) = Yl{y-rrii). 


2=1 


2=1 


Q{'t,y) = 


(1 - ytff 


life + t — mi + ytrfii), 


nnm{Q{t,y)) = JJ(2/(1 + tmi) + t - m*). 

2=1 

Since p{y) and num(Q(t, ?/)) have both positive degree in y, P{t) = 0 if and 
only if they have a common factor with positive degree in y (see Chapter 3 
in ID)- In turn, this happens if and only if there exists r G {1, ..., n} such 
that {y — mr) divides num((5(t, ?/)). As a consequence num(Q(t,mr)) = 0. 
Then there exists j G {1,... ,n} with (m^ + t — rfij + tmrfhj) = 0, and so 
mr = fhj, 1 + fhj mr = 0 hold. Since p{y) and q{y) are real polynomials, this 
implies that |/^ + 1 divides both polynomials. □ 

Finally let us consider the case cos 6 = 0. Here we have a = ip,, where 
p = \a\ ■ sin6*, and therefore a = —ip. By using b = ri{a,a) and plugging 
these relationships into the system S, we reach a polynomial system with just 
one variable, p. Hence, here we just need to check whether or not the gcd 
of the (univariate, and with complex coefficients) polynomials in the system 
has some nonzero root. 

Hence, we have the following algorithm Impl-Sim-General to check sim¬ 
ilarity in the general case. 
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Algorithm Impl-Sim-General 

Require: Two irreducible, implicit algebraic curves f{x,y), g{x,y), neither 
lines nor circles, satisfying the hypotheses of the general case. 

Ensure: Whether or not they are related by an orientation-preserving sim¬ 
ilarity. 

1: Find the system S. 

2: Find j such that ctn-j,j 7^ 0 and Aj 7^ 0 . 

3: Find A (a, a) from ([^. 

4: Find b = r]{a,a) b = f]{a,a) and from ([^ and Q. 

5: Substitute A := A(a,a), b := ri{a,a), b = fj{a,a) in the system S. 
[Similarities with cos 9 7^ 0:] . 

6: Compute the polynomial P{t). 

7 : Substitute a = r ■ (1 -|- iuj) , a = r ■ (1 — iuj) , b = ri{a,a), b = fi{a,a) 
in S. If P{t) is not identically zero, add the equation P{oj) = 0 to the 
system. 

8: Check if the system has any solution with r G M, r 7^ 0 and a; G M. 

9: In the affirmative case, return “The curves are related by an 

orientation-preserving similarity”. 

[Similarities with cos 6 = 0.] 

10: Substitute a = iy, a = —iy, b = ri{a,a), b = fi{a,a) in S. 

11: Check if the univariate polynomials in y obtained this way have a non¬ 
trivial gcd with some nonzero real root. 

12: In the affirmative case, return “The curves are related by an 

orientation-preserving similarity”. 

13: If no similarities have been found, return “The curves are not 
related by an orientation-preserving similarity”. 
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5. The special case. 

From Lemmain Section]^ here we have j = 0. Furthermore, an,o 7 ^ 0 
and f3nfl ^ 0. Now if a„_i,o ^ 0, then from Q with j = 0 we can write a 
as a linear function of b, b, i.e. a = p{b, b). Let us see that this requirement 
can always be fulhlled by means of an appropriate translation. 

Lemma 6 . Let T{z) = z + k represent a translation by k & C, and let 
Cl = T(Ci). Also, let j G { 0 ,... ,n — 1}. The coefficient of Ci of bidegree 
{n-j- 1 ,j) is 

(Xfi—j—lJ J —Ij 'L K • (n j K • (^j (12) 


As a consequence, if On-jj 7 ^ 0 and an-j-ij = 0 then for almost all k G C, 
we have otn-j-ij 7 ^ 0 . 


Proof. The coefficient of Ci = T[Ci) of bidegree (n — j — 1, j) comes from the 
terms of F with bidegrees {n — j — 1, j) and — j — + Each 

of these terms provides, in turn, the terms of (13): the term of F of bidegree 
[n — j,j) provides the hrst term of (13), the term of bidegree {n — j — l,j) 
provides the term in k, and the term of {n — j — l,j + 1) provides the term 
in K. As for the second part of the statement, assume that an-jj 7 ^ 0, 
an-j-i,j = 0, and suppose by contradiction that An-j-ij vanishes for all n. 
By substituting k = 1 in (13), we have (n — j)an-j,j + {j + l)Q;n-j-i,j+i = 0; 
by substituting k = i in (13), we get that {n—j)an-jj — {j + l)an-j-ij+i = 0. 
Putting these two equations together, we deduce that an-j-ij+i = ctn-j,j = 
0 , contradicting that otn-jy 7 ^ 0 . □ 


Since similarities form a group, Ci and C 2 are similar if and only if 
Cl = T{Ci) and C2 are similar. Therefore, there is no problem in apply¬ 
ing a translation T on Ci. Furthermore, by Lemma a random translation 
will suffice. Now after replacing Ci by T(Ci), we can write a = ^{b,b), where 
^ is linear. By substituting a = ^{b, b) into the system S, we are led to a bi¬ 
variate polynomial system S* with complex coefficients, where the unknowns 
are b, b. In turn, this gives rise to a real bivariate polynomial system where 
the unknowns are the real and imaginary parts 61,62 of b. Then the curves 
are similar if and only if this system has a real solution where a 7 ^ 0 . 

Hence, we have the following algorithm Impl-Sim-Spec to check similar¬ 
ity in this case. 
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Algorithm Impl-Sim-Spec 

Require: Two irreducible, implicit algebraic curves f{x,y), g{x,y), neither 
lines nor circles, satisfying the hypotheses of the special case. 

Ensure: Whether or not they are related by an orientation-preserving sim¬ 
ilarity. 

1: Check whether or not an-i,o 7^ 0. 

2: If Un-ifl = 0 then hnd a translation T{z) = z + k such that dn-i,i 7^ 0, 
and replace Ci —)■ T(Ci). 

3: Write a = .^(h, h), A (a, a). 

4: Substitute the above expressions into S. 

5: Write b = bi + ib 2 to get the real bivariate polynomial system S*. 

6: Check if S* has a real solution with a ^ 0. 

7: In the affirmative case, return “The curves are related by 
an orientation-preserving similarity”. Otherwise, return 

“The curves are not related by an orientation-preserving 
similarity”. 


6. Observations on the “orientation-reversing” case. 

In order to look for orientation-reversing similarities transforming Ci into 
C 2 , we must check the consistency of the system S' (see the beginning of 
Section]^, which stems from Equation ([^. The analysis is very similar to 
orientation-preserving similarities. However, for completeness we summarize 
here some observations that must be taken into account. 

First, Equation ([^ must be replaced by 



(14) 


which is deduced by comparing the terms of z^z'^ ^ at both sides of (|^. 
Similarly, Equation ([^ must be replaced by 


(n ■ h -|- (j T l)/3jj_j_i j_|_i ■ b j 

and (|^ must be replaced by 




(j + ■ h + (n - j)f3n-jj ■ b + /3n-j-l,j = 
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Lemma 1^ also holds here, and again we mnst distingnish between the gen¬ 
eral case and the special case, which are formnlated exactly in the same way. 
There are some slight changes in the discnssion at the beginning of Section]^ 
thongh. In this case, the curve defined by /„(x, ?/) = 0 is transformed into the 
curve defined by gn{x,y) = 0 under the transformation h{z) = az, which ge¬ 
ometrically is equivalent to a symmetry with respect to the x-axis, followed 
by a rotation about the origin by an angle 9. Again, the case cos( 6 *) = 0 
must be considered separately, and leads to a simpler case, which presents 
no new difficulties. The case (1) in Section leads to either P{t) = t, or 
P(t) = t{t — d/c) or P{t) = t — d/c. In case (2) , we have 6 * = 71 -|- 72 . So 
denoting m = tg( 7 i), fh = tg( 72 ), and since 72 = 6 * — 7 i, we get that: 

oj — m 

m = -. 

1 -|- mu 


Hence, a; is a root of 

P{t) = Resy |^p(|/),numer (^q = 0- 

Furthermore, Proposition]^ also holds here. Additionally, in the special case. 
Lemma 1^ is applied to 

7. Implementation and timings 

Algorithms Impl-Sim-General and Impl-Sim-Spec have been imple¬ 
mented in Maple 2015 m)- The code is available in [6]. Since we used 
Maple’s generic solver SolveTools:-PolynomialSystem, the algorithm is 
deterministic. 

7.1. Example I 

Let Ci,C 2 be defined by f{x,y) and g{x,y), where 

f{x, y) = 15 x^y — 40 xy"^ — Iby'^ + b x^ + b xy —‘iby'^ + bx — by + 2 


and 


g{x, y) = y^ + 2 xy'^ — x'^y — xy — 2x^ -\-1. 


Both curves are shown in Figure]^ Using Lemmaone can check that we 
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Figure 1: The curve Ci (left) and a related curve C 2 obtained by similarity 
(right). 


fall into the general case. By applying the algorithm Impl-Sim-General, we 
get an orientation-preserving similarity h(z) = az -h b with 

P(t) = -125 (t^+2t^-t-2) (448f®+4416t®+8880t^-1920f^-8880f2+4416f-448). 


Adding the equation P(uj) = 0 to the system in r, uj obtained from S, we get 


. . , 17r 19ujr 

a = Kl + «c.) = l-2,. 

The ratio of this similarity is equal to ^/5. 


1 Ticnr Hr 1, 

-h(-^-h-)i = 1- 

10 ^ 100 50 5^ 


■i. 


1.2. Example II 

Let Ci,C 2 be defined by f{x,y) and g{x,y), where 
f{x, y) = x'^ + 2 x^y"^ + y'^ — % x‘^y —‘&y^ + 12x‘^ — Q xy + 29y‘^ + 12x — IQy., 
and 

g{x,y) =2x'^ + Ax^y"^ + 2y^ — x^ + y^. 

Both curves are shown in Figure In this case we also fall into the gen- 
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Figure 2: A lemniscate Ci (left) and a related lemniscate C 2 obtained by 
similarity (right). 


eral case. However, here P(t) is identically zero. By applying the algo¬ 
rithm Impl-Sim-General, we hrst hnd two orientation-preserving similari¬ 
ties h(z) = az -h b. In terms of uj and r, we get b = —2ir (1 + ioj) and 
A = 2r^ (1 -|- iojY (1 — iojY. The hrst orientation-preserving similarity corre¬ 
sponds to a; = —3 and r = 1/10. Since a = r(l -|- ica), we have 


1 




b 



The second orientation-preserving similarity corresponds to uj 
r = —1/10, and 


—3 and 


1 3 . , 3 1 . 

a = - 1- i, 0 = —I— i, 

10 10 ’ 5 5’ 


1 


Additionally we hnd two orientation-reversing similarities h{z) = az + b. In 
terms of u and r, we have b = 2ir {1 + iu) and A = 2r^ (1 -|- iu)‘^ (1 — 

The hrst orientation-reversing similarity corresponds to a; = 3 and r = 1/10, 
and 



3 1 

■5 + 5'’ 


1 


The second orientation-reversing similarity corresponds to u 
— 1/10, and 



1 


3 and r = 
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Figure 3: The curve Ci (left) and a related curve C 2 obtained by similarity 
(right). 


The ratio of all these similarities is equal to 
7.5. Example III 

Let Ci,C 2 be defined by f{x,y) and g{x,y), where 
f{x,y) = 19x^+90x^y—18xy^+35y^+51x^+237xy—90y^+39x+195y—l, 
and 

g{x,y) = x^ + y^ — 3yx (Descartes’ Folium). 

Both curves are shown in Figure Using Lemma we observe that we 
fall into the special case addressed in Section By applying the algorithm 
Impl-Sim-Spec, we first find an orientation-preserving similarity h{z) = az+ 
b, where b = 3 — Ai, a = 3 — 2i, X = 1. Additionally, we hnd an orientation- 
reversing similarity h{z) = az + b, where b = —4 + 3i, a = —2 -|- 3z, A = 1. 
The ratio of these similarities is equal to -\/l3- 

7.4- Performance 

Algorithms Impl-Sim-General and Impl-Sim-Spec can be unified into 
one algorithm that first checks the case we are in, and then executes Al¬ 
gorithms Impl-Sim-General or Impl-Sim-Spec accordingly. We tested the 
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t 


r = 1 

r = 2 

r = 4 

r = 8 

r = 16 

r = 32 

d = 

5 

0.70 

0.73 

0.73 

0.79 

0.82 

2.73 

d = 

6 

1.04 

1.11 

1.22 

1.51 

5.37 

9.11 

d = 

7 

1.84 

3.24 

8.10 

11.10 

12.76 

14.71 

d = 

8 

15.94 

17.80 

20.40 

21.91 

25.23 

32.40 

d = 

9 

24.82 

26.86 

27.79 

30.75 

36.20 

42.15 

d = 

10 

35.84 

38.60 

41.81 

47.45 

52.13 

72.27 


Table 1: Average CPU time t (seconds) of the algorithms applied to random 
dense bivariate polynomials. 


performance of this algorithm on curves of various degrees d and bitsize 
More precisely, for every d G |5..10] and t = 2\i & [0..5], we generated 20 
random dense bivariate polynomials g{x, y) of degree d and bitsizes bounded 
by T. Once generated, we computed F{z, z) := G(az + b, az + b), with (a, b) 
a pair of random complex numbers with coefficients in the range |—10..10]. 
Finally we tested the algorithms with F(z, z) and G(z, z). Observe that even 
if the bitsize of g(x, y) is bounded by r, the bitsize of F{z, z) is usually bigger 
than r. The code, as well as the random inputs, are given in [B]; the timings 
are provided in Table [T} 

After examining the correlation matrix of the data we observed that the 
degree highly influences the timing. The influence of the bitsize is weaker. By 
using Maple’s command Statistics: - NonlinearFit, we compared several 
following multivariate model functions and found that the following function 
hts our data quite well, 

t = -0.0987387345679027195 d^ + 3.89759623015877255 
+0.180291478322672782 rd- 33.7316950021062496 d 
-0.853707190239283653 r + 84.1154055239685192. 


Figure l^shows the behavior of this model, and compares it with the data col- 


•^The bitsize t := |"log 2 fc] + 1 of an integer k is the number of bits needed to represent 
it. When we refer to the bitsize of an implicit algebraic curve, we imply the maximum of 
the bitsizes of the absolute values of the coefficients of the implicit equation. 
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lected in our experiments. Notice that this model suggests a time complexity 
0{d^ + dr) for the algorithm. 


° deg.=5 ° deg.=6 ° deg.=7 “ deg.=8 


° T=1 ° t=2 ° 1=4 “ t=8 

° deg.=9 ° deg.=10 


□ x = \6 “ t = 32 




Figure 4: Degrees, bitsizes and timings. 


We must emphasize that these timings corresponds to dense polynomi¬ 
als. Therefore, in practice the performance is better than Table [T] predicts. 
In order to illustrate this statement we tested the algorithms with several 
classical curves, which we list below: 

• Folium of Descartes, 

g{x,y) = -3yx + y^, 

• Lemniscate of Bernoulli, 

g{x,y) =2x'^ + 4 x‘^y‘^ + 2y‘^ — x^ + y'^, 


• Epitrochoid, 

g{x,y) = x^ + 2x‘^y‘^ -I- — 34— 34 -|- 96 a; — 63, 
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• Offset curve to the cardioid, 

g{x,y) = + Ax^y'^ + + 16x^y + 

+48 x^y^ + 16y’^ - 140 x^ - 324 x^y"^ 

-228 a; V - 44/ - 1552 x^y - 2848 a;^/ 

-1296/ + 2416 a;^ - 864 a; V - 3024/ + 28800 x‘^y 
+21888 / + 43776 a;^ + 124416 / + 228096 y + 145152, 


• Hypocycloid, 

g{x,y) = a;® + 4a;®|/^ + 6a;‘^/+ 4a;^l/® +/l2a;® + 36a;^l/^ + 36a;^/ 
+ + 12 / - 512 x^ + 5120 x^y^ - 2560 a;/ + 150 x^ 

+300 x^y"^ + 150 / + 67500 x^ + 67500 / - 759375, 


• 4-^eaf Rose, 

g{x, y) = x^ + 3 a;^/ + 3 a;^/ + / — a;^ + 2 x^y"^ — /, 

• 8-leaf Rose, 

g{x,y) = a;^° + 5 a;®/+ 10 a;®/+ lOa;"^/ + 5a;^/+ — a;® + 12a;®l/^ 

—38a;^/ + 12 a;^/ — /, 

• 12-leaf Rose, 

g{x,y) = a;^^ + 7a;^^/ + 21a;^°/ + 35a;®|/® + 35 a;®/ + 21a;^l/^° 

+7 + y^^ — x^‘^ + 30 x^^y'^ — 255 a;®/ + 452 a;®?/® 

-255a;V + 30a;2|/i0-/2, 

• 16-leaf Rose, 

g{x,y) = a;^® + 9a;^®l/^ + 36a;^^/+ 84a;^^?/® + 126a;^®l/® + 126a;®|/^® 
+84 a;®!/^^ + 36 + 9 x'^y^^ + y^® — a;^® + 56 a;^^/ 

-924a;^2/ + 3976 a;^®/ - 6470 a;®/ + 3976 a;®?/^® - 924 a ;^|/^2 
+56 — y^®. 


22 



Curve 

Descartes’ 

Bernoulli’s 

Epitrochoid 

Cardioid 

Hypocycloid 


folium 

lemniscate 


offset 




CX3 

G 

0 

G 

degree 

3 

4 

4 

8 

8 

CPU time 

0.46 

0.39 

0.46 

2.62 

2.37 

Curve 

4-leaf rose 

8-leaf rose 

12-leaf rose 

16-leaf rose 

20-leaf rose 





# 

# 

degree 

6 

10 

14 

18 

22 

CPU time 

0.75 

0.79 

1.56 

5.67 

9.85 


Table 2: Average CPU time (seconds) of the algorithms for well-known 
curves. 


• 20-leaf Rose, 

g{x,y) = -t--f -I-330x^^l/®-I- 

-1-462 -I- 330 165 x^y^^ -\- 55 x'^y^^ -\-11 -|- 

-x2° + 90x^V - 2445 + 19320 x^V - 63090x^2/ 

+92252 x^°?/^° - 63090 x®2/^^ + 19320 - 2445 

+90xV® 

The timings corresponding to these curves are given in Table 
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